Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans.

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Tungsten-anisole complex provides 3,6-substituted cyclohexenes for highly diversified chemical libraries.

Medicinal chemists use vast combinatorial molecular libraries to develop leads for new pharmaceuticals. The syntheses of these compounds typically rely on coupling molecular fragments through atoms with planar (sp2) geometry. These so-called flat molecules often lack the protein binding site specificity needed to be an effective drug. Here, we demonstrate a coupling strategy in which a cyclohexene is used as a linker to connect two diverse molecular fragments while forming two new tetrahedral (sp3) stereocenters. These connections are made with the aid of a tungsten complex that activates anisole toward an unusual double protonation, followed by sequential nucleophilic additions. As a result, either cis- or trans-disubstituted cyclohexenes can be prepared with a range of chemical diversity unparalleled by other dearomatization methods.

Integrative single-cell characterization of a frugivorous and an insectivorous bat kidney and pancreas.

Frugivory evolved multiple times in mammals, including bats. However, the cellular and molecular components driving it remain largely unknown. Here, we use integrative single-cell sequencing (scRNA-seq and scATAC-seq) on insectivorous (Eptesicus fuscus; big brown bat) and frugivorous (Artibeus jamaicensis; Jamaican fruit bat) bat kidneys and pancreases and identify key cell population, gene expression and regulatory differences associated with the Jamaican fruit bat that also relate to human disease, particularly diabetes. We find a decrease in loop of Henle and an increase in collecting duct cells, and differentially active genes and regulatory elements involved in fluid and electrolyte balance in the Jamaican fruit bat kidney. The Jamaican fruit bat pancreas shows an increase in endocrine and a decrease in exocrine cells, and differences in genes and regulatory elements involved in insulin regulation. We also find that these frugivorous bats share several molecular characteristics with human diabetes. Combined, our work provides insights from a frugivorous mammal that could be leveraged for therapeutic purposes.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans.

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using ATAC-STARR-seq. In addition to thousands of cis changes, we discover an unexpected number (~10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie >50% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-to regulatory differences between species.

Function and Constraint in Enhancer Sequences with Multiple Evolutionary Origins.

Thousands of human gene regulatory enhancers are composed of sequences with multiple evolutionary origins. These evolutionarily "complex" enhancers consist of older "core" sequences and younger "derived" sequences. However, the functional relationship between the sequences of different evolutionary origins within complex enhancers is poorly understood. We evaluated the function, selective pressures, and sequence variation across core and derived components of human complex enhancers. We find that both components are older than expected from the genomic background, and complex enhancers are enriched for core and derived sequences of similar evolutionary ages. Both components show strong evidence of biochemical activity in massively parallel report assays. However, core and derived sequences have distinct transcription factor (TF)-binding preferences that are largely similar across evolutionary origins. As expected, given these signatures of function, both core and derived sequences have substantial evidence of purifying selection. Nonetheless, derived sequences exhibit weaker purifying selection than adjacent cores. Derived sequences also tolerate more common genetic variation and are enriched compared with cores for expression quantitative trait loci associated with gene expression variability in human populations. In conclusion, both core and derived sequences have strong evidence of gene regulatory function, but derived sequences have distinct constraint profiles, TF-binding preferences, and tolerance to variation compared with cores. We propose that the step-wise integration of younger derived with older core sequences has generated regulatory substrates with robust activity and the potential for functional variation. Our analyses demonstrate that synthesizing study of enhancer evolution and function can aid interpretation of regulatory sequence activity and functional variation across human populations.

Modulated Calcium Homeostasis and Release Events Under Atrial Fibrillation and Its Risk Factors: A Meta-Analysis.

Background: Atrial fibrillation (AF) is associated with calcium (Ca2+) handling remodeling and increased spontaneous calcium release events (SCaEs). Nevertheless, its exact mechanism remains unclear, resulting in suboptimal primary and secondary preventative strategies. Methods: We searched the PubMed database for studies that investigated the relationship between SCaEs and AF and/or its risk factors. Meta-analysis was used to examine the Ca2+ mechanisms involved in the primary and secondary AF preventative groups. Results: We included a total of 74 studies, out of the identified 446 publications from inception (1982) until March 31, 2020. Forty-five were primary and 29 were secondary prevention studies for AF. The main Ca2+ release events, calcium transient (standardized mean difference (SMD) = 0.49; I 2 = 35%; confidence interval (CI) = 0.33-0.66; p < 0.0001), and spark amplitude (SMD = 0.48; I 2 = 0%; CI = -0.98-1.93; p = 0.054) were enhanced in the primary diseased group, while calcium transient frequency was increased in the secondary group. Calcium spark frequency was elevated in both the primary diseased and secondary AF groups. One of the key cardiac currents, the L-type calcium current (ICaL) was significantly downregulated in primary diseased (SMD = -1.07; I 2 = 88%; CI = -1.94 to -0.20; p < 0.0001) and secondary AF groups (SMD = -1.28; I 2 = 91%; CI = -2.04 to -0.52; p < 0.0001). Furthermore, the sodium-calcium exchanger (INCX) and NCX1 protein expression were significantly enhanced in the primary diseased group, while only NCX1 protein expression was shown to increase in the secondary AF studies. The phosphorylation of the ryanodine receptor at S2808 (pRyR-S2808) was significantly elevated in both the primary and secondary groups. It was increased in the primary diseased and proarrhythmic subgroups (SMD = 0.95; I 2 = 64%; CI = 0.12-1.79; p = 0.074) and secondary AF group (SMD = 0.66; I 2 = 63%; CI = 0.01-1.31; p < 0.0001). Sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) expression was elevated in the primary diseased and proarrhythmic drug subgroups but substantially reduced in the secondary paroxysmal AF subgroup. Conclusions: Our study identified that ICaL is reduced in both the primary and secondary diseased groups. Furthermore, pRyR-S2808 and NCX1 protein expression are enhanced. The remodeling leads to elevated Ca2+ functional activities, such as increased frequencies or amplitude of Ca2+ spark and Ca2+ transient. The main difference identified between the primary and secondary diseased groups is SERCA expression, which is elevated in the primary diseased group and substantially reduced in the secondary paroxysmal AF subgroup. We believe our study will add new evidence to AF mechanisms and treatment targets.

Modeling the Evolutionary Architectures of Transcribed Human Enhancer Sequences Reveals Distinct Origins, Functions, and Associations with Human Trait Variation.

Despite the importance of gene regulatory enhancers in human biology and evolution, we lack a comprehensive model of enhancer evolution and function. This substantially limits our understanding of the genetic basis of species divergence and our ability to interpret the effects of noncoding variants on human traits. To explore enhancer sequence evolution and its relationship to regulatory function, we traced the evolutionary origins of transcribed human enhancer sequences with activity across diverse tissues and cellular contexts from the FANTOM5 consortium. The transcribed enhancers are enriched for sequences of a single evolutionary age ("simple" evolutionary architectures) compared with enhancers that are composites of sequences of multiple evolutionary ages ("complex" evolutionary architectures), likely indicating constraint against genomic rearrangements. Complex enhancers are older, more pleiotropic, and more active across species than simple enhancers. Genetic variants within complex enhancers are also less likely to associate with human traits and biochemical activity. Transposable-element-derived sequences (TEDS) have made diverse contributions to enhancers of both architectures; the majority of TEDS are found in enhancers with simple architectures, while a minority have remodeled older sequences to create complex architectures. Finally, we compare the evolutionary architectures of transcribed enhancers with histone-mark-defined enhancers. Our results reveal that most human transcribed enhancers are ancient sequences of a single age, and thus the evolution of most human enhancers was not driven by increases in evolutionary complexity over time. Our analyses further suggest that considering enhancer evolutionary histories provides context that can aid interpretation of the effects of variants on enhancer function. Based on these results, we propose a framework for analyzing enhancer evolutionary architecture.

Genetic risk for major depressive disorder and loneliness in sex-specific associations with coronary artery disease.

Major depressive disorder (MDD) and loneliness are phenotypically and genetically correlated with coronary artery disease (CAD), but whether these associations are explained by pleiotropic genetic variants or shared comorbidities is unclear. To tease apart these scenarios, we first assessed the medical morbidity pattern associated with genetic risk factors for MDD and loneliness by conducting a phenome-wide association study in 18,385 European-ancestry individuals in the Vanderbilt University Medical Center biobank, BioVU. Polygenic scores for MDD and loneliness were developed for each person using previously published meta-GWAS summary statistics, and were tested for association with 882 clinical diagnoses ascertained via billing codes in electronic health records. We discovered strong associations with heart disease diagnoses, and next embarked on targeted analyses of CAD in 3893 cases and 4197 controls. We found odds ratios of 1.11 (95% CI, 1.04-1.18; P 8.43 × 10-4) and 1.13 (95% CI, 1.07-1.20; P 4.51 × 10-6) per 1-SD increase in the polygenic scores for MDD and loneliness, respectively. Results were similar in patients without psychiatric symptoms, and the increased risk persisted in females even after adjusting for multiple conventional risk factors and a polygenic score for CAD. In a final sensitivity analysis, we statistically adjusted for the genetic correlation between MDD and loneliness and re-computed polygenic scores. The polygenic score unique to loneliness remained associated with CAD (OR 1.09, 95% CI 1.03-1.15; P 0.002), while the polygenic score unique to MDD did not (OR 1.00, 95% CI 0.95-1.06; P 0.97). Our replication sample was the Atherosclerosis Risk in Communities (ARIC) cohort of 7197 European-ancestry participants (1598 incident CAD cases). In ARIC, polygenic scores for MDD and loneliness were associated with hazard ratios of 1.07 (95% CI, 0.99-1.14; P = 0.07) and 1.07 (1.01-1.15; P = 0.03), respectively, and we replicated findings from the BioVU sensitivity analyses. We conclude that genetic risk factors for MDD and loneliness act pleiotropically to increase CAD risk in females.

Accounting for diverse evolutionary forces reveals mosaic patterns of selection on human preterm birth loci.

Currently, there is no comprehensive framework to evaluate the evolutionary forces acting on genomic regions associated with human complex traits and contextualize the relationship between evolution and molecular function. Here, we develop an approach to test for signatures of diverse evolutionary forces on trait-associated genomic regions. We apply our method to regions associated with spontaneous preterm birth (sPTB), a complex disorder of global health concern. We find that sPTB-associated regions harbor diverse evolutionary signatures including conservation, excess population differentiation, accelerated evolution, and balanced polymorphism. Furthermore, we integrate evolutionary context with molecular evidence to hypothesize how these regions contribute to sPTB risk. Finally, we observe enrichment in signatures of diverse evolutionary forces in sPTB-associated regions compared to genomic background. By quantifying multiple evolutionary forces acting on sPTB-associated regions, our approach improves understanding of both functional roles and the mosaic of evolutionary forces acting on loci. Our work provides a blueprint for investigating evolutionary pressures on complex traits.

Chronic interleukin-1 exposure drives haematopoietic stem cells towards precocious myeloid differentiation at the expense of self-renewal.

Haematopoietic stem cells (HSCs) maintain lifelong blood production and increase blood cell numbers in response to chronic and acute injury. However, the mechanism(s) by which inflammatory insults are communicated to HSCs and their consequences for HSC activity remain largely unknown. Here, we demonstrate that interleukin-1 (IL-1), which functions as a key pro-inflammatory 'emergency' signal, directly accelerates cell division and myeloid differentiation of HSCs through precocious activation of a PU.1-dependent gene program. Although this effect is essential for rapid myeloid recovery following acute injury to the bone marrow, chronic IL-1 exposure restricts HSC lineage output, severely erodes HSC self-renewal capacity, and primes IL-1-exposed HSCs to fail massive replicative challenges such as transplantation. Importantly, these damaging effects are transient and fully reversible on IL-1 withdrawal. Our results identify a critical regulatory circuit that tailors HSC responses to acute needs, and is likely to underlie deregulated blood homeostasis in chronic inflammation conditions.

Anti-EFNA4 Calicheamicin Conjugates Effectively Target Triple-Negative Breast and Ovarian Tumor-Initiating Cells to Result in Sustained Tumor Regressions.

PURPOSE: Triple-negative breast cancer (TNBC) and ovarian cancer each comprise heterogeneous tumors, for which current therapies have little clinical benefit. Novel therapies that target and eradicate tumor-initiating cells (TIC) are needed to significantly improve survival. EXPERIMENTAL DESIGN: A panel of well-annotated patient-derived xenografts (PDX) was established, and surface markers that enriched for TIC in specific tumor subtypes were empirically determined. The TICs were queried for overexpressed antigens, one of which was selected to be the target of an antibody-drug conjugate (ADC). The efficacy of the ADC was evaluated in 15 PDX models to generate hypotheses for patient stratification. RESULTS: We herein identified E-cadherin (CD324) as a surface antigen able to reproducibly enrich for TIC in well-annotated, low-passage TNBC and ovarian cancer PDXs. Gene expression analysis of TIC led to the identification of Ephrin-A4 (EFNA4) as a prospective therapeutic target. An ADC comprising a humanized anti-EFNA4 monoclonal antibody conjugated to the DNA-damaging agent calicheamicin achieved sustained tumor regressions in both TNBC and ovarian cancer PDX in vivo. Non-claudin low TNBC tumors exhibited higher expression and more robust responses than other breast cancer subtypes, suggesting a specific translational application for tumor subclassification. CONCLUSIONS: These findings demonstrate the potential of PF-06647263 (anti-EFNA4-ADC) as a first-in-class compound designed to eradicate TIC. The use of well-annotated PDX for drug discovery enabled the identification of a novel TIC target, pharmacologic evaluation of the compound, and translational studies to inform clinical development.

Re-entry into quiescence protects hematopoietic stem cells from the killing effect of chronic exposure to type I interferons.

Type I interferons (IFN-1s) are antiviral cytokines that suppress blood production while paradoxically inducing hematopoietic stem cell (HSC) proliferation. Here, we clarify the relationship between the proliferative and suppressive effects of IFN-1s on HSC function during acute and chronic IFN-1 exposure. We show that IFN-1-driven HSC proliferation is a transient event resulting from a brief relaxation of quiescence-enforcing mechanisms in response to acute IFN-1 exposure, which occurs exclusively in vivo. We find that this proliferative burst fails to exhaust the HSC pool, which rapidly returns to quiescence in response to chronic IFN-1 exposure. Moreover, we demonstrate that IFN-1-exposed HSCs with reestablished quiescence are largely protected from the killing effects of IFNs unless forced back into the cell cycle due to culture, transplantation, or myeloablative treatment, at which point they activate a p53-dependent proapoptotic gene program. Collectively, our results demonstrate that quiescence acts as a safeguard mechanism to ensure survival of the HSC pool during chronic IFN-1 exposure. We show that IFN-1s can poise HSCs for apoptosis but induce direct cell killing only upon active proliferation, thereby establishing a mechanism for the suppressive effects of IFN-1s on HSC function.